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1.0 


INTRODUCTION 


The  ability  to  discriminate  at  long  range  between  sound 
sources  that  are  at  different  ocean  depths  would  be  useful  for  a 
variety  of  purposes.  Ocean-acoustic  experiments,  underwater 
communication,  and  acoustic  detection  of  submarines  often  suffer 
from  background  noise  originating  at  the  ocean  surface,  either  from 
surface  waves  or  from  surface  ships.  If  an  acoustic  source,  whose 
purpose  may  be  either  ocean  monitoring  or  cummunication,  can  be 
located  several  hundred  meters  below  the  surface,  then  a  vertical 
array  might  be  used  to  discriminate  against  unwanted  noise. 

A  naive  investigator  might  argue  against  the  possibility  of 
such  vertical  discrimination.  The  argument  goes  as  follows:  We  are 
interested  only  in  long-range  systems  (say  500  km  as  an  example), 
which  restricts  us  to  low  frequencies  (say  150  Hz  as  an  example) 
because  of  absorption.  A  simple  calculation  of  the  resolution  of  an 
array  would  say  that  to  distinguish  two  sources  separated  by  50  m  at 
1000  km  range  would  require  a  200-km-long  vertical  array.  Since  the 
depth  of  the  ocean  is  less  than  5  km,  the  task  is  Impossible. 

Yet,  several  experiments  by  Fisher  and  Williams using  a 
vertical  array  at  the  depth  of  the  sound  channel  axis  have  indicated 
a  substantial  difference  in  vertical  arrival  angle  of  the  energy 


arriving  from  two  sources  at  depths  of  25  and  80  meters.  These 
differences  persist  out  to  ranges  of  1500  kilometers  and  possibly 
more. 

Although  simple  ray-tracing  arguments  support  the 
experimental  results,  the  validity  of  ray  tracing  at  150  Hz  to 
ranges  greater  than  1000  km  is  open  to  question. 

It  is  the  purpose  of  this  note  to  present  full 
wave-propagation  calculations,  including  diffractive  effects  in  the 
parabolic  approximation,  that  quantitatively  predict  the  behavior  of 
finite  arrays,  even  near  caustics.  The  results  confirm  the 
experimental  observations,  and  provide  a  simple  way  to  predict 
results  for  a  variety  of  sound  channels.  The  calculations  were 
carried  out  with  parabolic-equation  numerical  simulation. 

A  simple  approximate  rule  for  vertical  discrimination  can 
be  given  as  follows:  The  approximate  length,  D  ,  of  a  vertical 
array  necessary  to  distinguish  between  two  acoustic  sources  is 

D  _  600  m  ^150  Hz  ^  ^ 

AT  f 

where  f  is  the  acoustic  frequency  and  AT  is  the  temperature 
difference  between  the  two  source  depths.  For  example,  the  CONTRK 


experiment  had  a  temperature  difference  of  about  4  °C  between  their 
25  m  and  80  m  source  depths,  so  would  only  have  required  a  150-m 
array  to  resolve  at  150  Hz.  Since  they  had  a  500  m  array,  they  did 
well.  The  important  point  is  that  (1)  shows  no  dependence  on 
range.  It  doesn't  matter  how  far  away  the  array  is  from  the  source; 
the  vertical  discrimination  is  undiminished.  Thus  the  presence  of 
the  sound  channel  completely  invalidates  the  naive  array-resolution 
argument  at  long  range. 

The  above  rule  applies  to  vertical  arrays  on  the 
sound-channel  axis.  Calculations  indicate  that  putting  an  array  at 
the  reciprocal  depth  of  one  of  the  sources  gives  very  poor 
discrimination  against  a  source  at  another  depth. 


2.0  EXPERIMENT  DESCRIPTION 

Both  in  1977  and  1979  a  vertical  array  was  hung  from  FLIP  a 
few  hundred  miles  off  San  Diego,  California.  The  configuration  is 
shown  if  Figure  1.  The  acoustic  sources  were  hung  from  the  R/V 
NEW  HORIZON,  which  steamed  at  about  six  knots  from  FLIP  is  a 
northwesterly  direction  (Figure  2). 

The  raytrace  for  acoustic  energy  from  the  80-m  source  is 
shown  in  Figure  3,  illustrating  the  dependence  of  the  peak  arrival 
angle  as  a  function  of  range,  including  shadow  zones  which  occur  out 
to  approximately  zone  5.  Some  experimental  results  that  verify  the 
features  of  Figure  3  are  shown  in  Figure  4.  More  detailed  angular 
distributions  as  a  function  of  range  are  shown  in  Figure  5.  The 
energy  arriving  at  the  array  is  shown  on  a  linear  scale. 

Plots  of  the  angular  distribution  averaged  over  range  at 
long  range  are  shown  in  Figures  6  and  7  for  the  25  m  and  80  m  deep 
sources,  respectively.  It  is  the  difference  between  Figures  6  and  7 
that  we  would  like  to  confirm  theoretically. 

Finally,  it  will  be  useful  to  have  the  angular  distribution 
of  ambient  noise  in  the  ocean  at  the  time  of  this  experiment. 
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3.0  PARABOLIC-EQUATION  PRINCIPLES 

It  is  known  that  at  150  Hz  the  propagation  of  sound  in  the 
ocean  volume  is  controlled  by  the  wave-speed  (index-of-ref raction) 
variations  that  occur  due  to  temperature,  salinity,  or  pressure 
difference  in  the  water.  The  most  pronounced  of  these  effects  is 
the  sound  channel  itself,  which  has  vertical  scales  of  variation  of 
hundreds  of  meters.  If  the  scales  of  variation  are  large  compared 
with  a  wavelength,  then  the  sound  is  deviated  from  its  original 
direction  by  only  a  small  angle.  Energy  within  the  sound  channel 
occupies  a  vertical-angle  region  within  +15°  from  the  horizontal; 
other  angles  are  absorbed  in  the  ocean  bottom. 

Under  these  conditions  of  small-angle  deviation,  the  sound 
propagation  obeys  a  parabolic  approximation  to  the  full  wave 
equation.  This  approximation  is  of  the  form 

"  lq  3zz*  +  *  i3xf  (2) 

where  f(x,z)  is  the  demodulated  sound  amplitude;  x  is  the 
range;  z  is  the  depth,  q  is  the  acoustic  wavenumber,  and  p(z) 
is  the  fractional  deviation  of  sound  speed  from  the  reference  value 


used  to  calculate  q  . 


The  boundary  condition  at  the  ocean  surface  is  a 


pressure-release  reflection 

f(x,0)  =0  (3 

and  the  boundary  condition  at  the  ocean  bottom  is  that  f(x,  -h)  b< 
only  an  outgoing  wave;  that  is,  a  perfectly  absorbing  ocean  bottom 
at  depth  h  . 

The  initial  wave  at  x  *  0  can  be  adjusted  to  give  an 
approximation  to  whatever  the  source  may  be:  for  example,  a  point 
source . 


The  computer  realization  of  the  PE  method  in  our  code  is  a 
split-step  Fast  Fourier  Transform. 


fCx+dx.z)  =  exp[y  iAdx]  exp[iBdx]  exp[y  iAdx]  f(x,z)  (V 


A  *  ~  3 
2q  zz 


(5; 


B  =  -qy(z) 


(6] 


and  the  evaluation  of  the  exponential  of  the  second  derivative  is 


carried  out  by  FFT: 


exp[iAdx]$(z)  =  F  ^expfiXdx]  F[(j»(z)]  } 


(7) 


A  =  -  1.  k2  (8) 

2q 

where  F  stands  for  Fourier  transform  from  z-space  to  k-space. 

The  surface  boundary  condition  is  treated  by  doubling  the 
depth  region  and  introducing  a  reflected  source  ¥(0,-z)  above  the 
ocean  surface. 

The  bottom  is  treated  by  multiplying  each  ^(x.z)  by  a 
depth-dependent  absorption  factor 

L(z)  =  exp  {-a  dx  exp[  (^-)2]  }  (9) 

where  our  runs  have  used  a  *  0.05  m  *  ,  8  =  0.04  h  ,  and 
1  *  5120  m.  This  allows  energy  above  a  depth  of  a  bit  more  than 
4  km  to  be  unaffected  by  the  bottom,  yet  gives  very  little 
reflection  from  the  bottom. 

The  FFT  is  carried  out  with  512  points  over  the  "real" 
ocean;  that  is  a  total  1024-point  transform  when  the  reflected  ocean 


is  included.  The  step  size  in  range  at  150  Hz  has  been  taken  as 
0.5  km. 


The  initial  wavef unction  ¥(0,z)  is 

z_2  Z+Z 

f(0,z)  -  0.9144(irX)“1/2{exp[-( - -)2  ]  -  exp[-( - -)  ]  }  (10) 

X  X 

where  X  is  the  acoustic  wavelength.  This  form  closely 
approximates  a  point  source  at  depth  zg  with  unit  intensity  at  a 
range  of  one  yard. 

The  intensity  of  the  wave  at  range  x  is  then 

I(x,z)  -  1 1  (11) 
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4.0 


PARABOLIC-EQUATION  SIMULATION  RESULTS 


The  detailed  sound-speed  profile  used  in  the  PE  simulations 
was  obtained  from  Bruce  Williams  as  representative  of  the  1979 
CONTRK  V  experiment  (Figure  9).  Three  source  depths  were  chosen: 

25  m,  80  m,  and  150  m.  It  is  seen  in  Figure  9  that  the  reciprocal 
depths  for  these  sources  are  3800  m,  2800  m,  and  1900  m, 
respectively. 

The  computer-simulated  data  were  analyzed  in  terms  of 
500-m-long  arrays  placed  at  two  depths;  the  first,  from  500-1000  m, 
simulated  the  CONTRK  V  experiment;  the  second,  from  3500-4000  m, 
simulated  a  reciprocal-depth  experiment.  The  computer  was  run  out 
to  a  range  of  1300  km,  and  data  was  kept  from  10-300  km  and  from 
1000-1300  km  in  one  kilometer  steps. 

4.1  Intensity  as  a  Function  of  Range  and  Depth 

Figures  10-12  show  |y|^  on  a  linear  scale  as  a  function 
of  range  and  depth.  The  first  convergence  zone  at  50-60  km  is  seen 
clearly,  as  is  the  absorbing  bottom.  The  traces  of  convergence 
zones  at  1000  km  are  seen  for  the  150-m-deep  source,  but  not  for  the 
25-m-deep  source.  These  figures  gives  confidence  that  the  code  is 
functioning  properly. 
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Intensity  as  a  Function  of  Range  and  Angle  with  an  Axis  Array 

2 

Figures  13-18  show  |y|  on  a  linear  scale  as  a  function 
of  range  and  vertical  angle,  where  zero  degrees  is  the  horizontal 
direction.  The  array  is  at  500-1000  m  depth,  straddling  the  sound 
axis . 

At  short  range  the  convergence  zones  are  seen  clearly,  and 
the  difference  between  the  angles  received  from  the  different 
sources  are  also  seen  clearly,  with  the  shallower  sources  giving 
steeper  angles. 

At  long  range  the  angle  differences  persist,  even  when  the 
vestiges  of  convergence  zones  completely  disappear.  It  is 
noteworthy  that  the  deeper  (150  m)  source  still  has  clearly  defined 
convergence  zones  even  at  1200  km  range.  The  steep  thermocline 
breaks  up  long-range  convergence  zones  for  sources  within  it,  but 
sources  below  it  will  still  have  strong  convergence  zones  beyond 
1000  km. 

The  angular  distributions,  averaged  over 
two-hundred-kilometer  range  increments,  are  shown  in  Figures 
19-24.  No  significant  differences  between  the  distributions  at 
0-200  km  and  those  at  1000-1200  km  are  seen,  but  the  difference  in 
angles  for  sources  at  different  depths  is  striking,  These  plots  can 


be  compared  with  the  experimental  data  in  Figures  6  and  7.  The 
comparison  is  quite  good. 

The  distribution  of  intensity  versus  range,  averaged  over 
all  vertical  angles,  but  using  the  axis  array,  is  shown  in  Figures 
25-30.  At  short  range  the  convergence  zones  are  so  marked  that 
detection  within  the  shadow  zones  is  likely  to  be  difficult.  At 
long  range,  the  energy  from  relatively  shallow  sources  is  spread 
more  or  less  uniformly,  so  that  gaps  in  detection  will  not  occur. 
However  the  150-m  source  still  shows  strong  convergence  zones  at 
1200  km,  as  remarked  previously. 

4.3  Intensity  as  a  Function  of  Range  and  Angle  with  a 

Reciprocal-Depth  Array 

It  has  been  suggested  that  an  array  at  the  reciprocal  depth 
for  a  given  source  may  have  good  discrimination  against  other 
sources.  The  reciprocal  depth  is  that  depth  that  has  the  same  sound 
speed  as  the  source  depth.  Thus  a  ray  leaving  the  source  with  zero 
angle  will  approach  the  reciprocal  depth  at  zero  angle. 

Figures  31-36  show  the  intensity  as  a  function  of  range  and 
vertical  angle  inpinging  on  an  array  at  3500-4000  depth,  which  is 
approximately  the  reciprocal  depth  of  a  25  m  source.  Convergence 
zones  at  short  range  are  seen,  and  differences  in  arrival-angle 
distributions  for  different  sources  are  seen,  but  the  fundamental 
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fact  is  that  the  energy  from  all  the  shallow  depths  is  concentrated 
o 

within  5  of  zero  angle,  so  that  discrimination  of  source  depth  is 
practically  nonexistent.  Without  a  doubt  the  axis  array  has  better 
discrimination. 

The  distributions  in  vertical  angle  averaged  over  range  are 
shown  in  Figures  37-42  and  the  distributions  in  range  averaged  over 
vertical  angle  are  shown  in  Figures  43-48. 
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CONCLUSIONS 


The  use  of  a  vertical  array  on  the  axis  to  distinguish  the 
depth  of  an  acoustic  source  at  long-range  and  low  frequency  is 
technically  feasible.  The  discrimination  power  of  the  array  depends 
on  the  shape  of  the  sound-channel,  and  is  independent  of  range  to  a 
good  approximation.  Parabolic-equation  computer  simulations  of  wave 
propagation  at  150  Hz  and  1300  km  have  confirmed  an  experimental 
result  on  this  topic. 

The  full-wave  propagation  results  indicate  that  a  simple 
rule  taken  from  ray  tracing  should  be  reasonably  accurate;  consider 
the  energy  emitted  by  an  acoustic  source  to  be  primarily  at  zero 
vertical  angle.  Then  the  vertical  angle  at  the  sound  channel  axis 
can  be  calculated  from  the  sound  speed  at  the  source  depth  and  the 
sound  speed  at  the  axis,  by  Snell's  Law.  The  difference  between  the 
axis  angles  of  sources  at  two  different  depths  is  then 

02  "  df  "  (  c(z2)  -  C(ZL)  )  02) 

o 

o 

In  the  usual  sound  channel  both  0^  and  0^  are  abour  12  ,  and  the 
difference  C(z^)  ~  C(z  )  is  dominated  by  the  temperature 
difference  AT  between  the  two  depths.  We  have 


C  1  AC  -  a  AT 
o 


(13) 


where  a  =  3  x  10-^  °C-^  .  The  resolution  of  a  vertical  array 
length  D  is  A/D  ,  so  that  the  array  length  D  needed  to  resolve 
two  sources  can  be  derived: 


D 


A9 

aAT 


(14) 


For  150  Hz  we  have  A  -  10  m  .  Taking  9  •  0.2  we  have 


D 


600 

AT 


(15) 


Obviously  two  sources  within  the  mixed  layer  will  not  be 
well  discriminated.  However  any  source  in  water  a  few  degrees 
cooler  than  surface  water  should  be  well  discriminated  by  a 
several-hundred-meter  array. 


A  study  of  the  depth  at  which  the  difference  in  sound-axis 

o 

vertical  angle  with  the  surface-limited  ray  is  equal  to  1  has  been 

done  by  J.  Lovett^.  Figures  49  and  50  show  the  depth  (in  feet)  at 

which  this  discrimination  would  be  possible  if  the  array  were  able 
o 

to  resolve  a  1  difference.  The  second  number  in  each  ocean  region 
is  the  axis  angle  of  the  surface-limited  ray. 
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As  expected,  those  regions  of  the  ocean  with  a  sharp 
thremocline  are  able  to  discriminate  smaller  depth  differences. 

The  use  of  a  deep  array  (either  moored  or  towed)  does  not 
look  attractive  for  discrimination  because  all  sources  put  a  broad 
range  of  angles  around  zero  degrees  into  such  an  array. 

The  signal-to-noise  gain  at  150  Hz  with  a  500  m  vertical 
array  can  be  calculated  by  comparing  the  width  of  the  angular 

o 

distribution  of  received  intensity,  which  appears  to  be  about  2  , 
with  the  «  30°  width  for  ambient  noise.  Thus  the  gain  should  be 
about  12  dB.  It  would  be  desirable  to  calculate  this  gain  with 
various  sound  channels,  and,  at  higher  acoustic  frequencies,  with 
account  taken  of  angular  spread  due  to  scattering  from  ocean 
fluctuations. 

One  difficulty  with  this  axis  array  is  that  shadow  zones 
are  significant,  either  for  short  ranges  or  deep  sources.  A  broad 
study  of  many  realistic  sound  channels  should  be  done  to  ascertain 
at  what  range  the  shadow  zones  fill  in.  This  will  depend  on  source 
depth.  Also  it  would  be  desirable  to  develop  an  analytic  approach 
beyond  raytracing  for  this  problem. 
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Figure  1.  Configuration  of  CONTRK  vertical  array. 
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Figure  2.  Path  of  New  Horizon  towing  the  acoustic  source  for  CONTRK  V. 


Figure  3.  Raytrace  for  acoustic  energy  from  the  80-m  source 
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Figure  5.  Received  intensity  as  a  function  of  vertical  angle  and  elapsed  time,  detected  by 
the  vertical  array. 
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Figure  16.  Information  from  a  dense  array  between  500  and  1000  m  depth. 
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Source:  150  m;  Long  Range:  Shallow  Receiver 
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Figure  24.  Information  from  a  dense  array  between  500  and  1000  m  depth 


Source:  25  m;  Short  Range:  Shallow  Receiver 
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Figure  25.  Information  from  a  dense  array  between  500  and  1000  m  depth. 
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Figure  29.  Information  from  a  dense  array  between  500  and  1000  m  depth. 
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Figure  36.  Information  from  a  dense  array  between  3500  and  4000  m  depth. 
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Figure  37.  Information  from  a  dense  array  between  3500  and  4000  m  depth. 
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I  Figure  39.  Information  from  a  dense  array  between  3500  and  4000  m  depth. 
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Figure  40.  Information  from  a  dense  array  between  3500  and  4000  m  depth. 
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Figure  42.  Information  from  a  dense  array  between  3500  and  4000  m  depth. 
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Source:  80  m;  Short  Range:  Deep  Receiver 

Figure  44.  Information  from  a  dense  array  between  3500  and  4000  m  depth. 
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Source:  25  m;  Long  Range;  Deep  Receiver 

Figure  46.  information  from  a  dense  array  between  3500  and  4000  m  depth 


Figure  49  Ocean  regions  The  first  number  is  the  depth  in  feet  at  which  a  horizontal  ray 
arrives  at  the  sound  channel  axis  at  1°  different  angle  than  the  surface  limited 
ray  The  second  number  is  the  axis  angle  of  the  surface-limited  array.  (See 
Reference  3.) 
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Figure  50.  Ocean  regions.  The  first  number  is  the  depth  in  feet  at  which  a  horizontal  ray 
arrives  at  the  sound  channel  axis  at  1°  different  angle  than  the  surface  limited 
ray.  The  second  number  is  the  axis  angle  of  the  surface-limited  array.  (See 
Reference  3.) 
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